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ABSTRACT 

We present near-infrared spectroscopy and narrow-band imaging at the wavelength of redshifted Ha for a 
sample of 30 high-redshift, far-infrared luminous galaxies. This sample is selected from surveys in the sub- 
millimeter, millimeter and radio wavebands and has complete redshift coverage with a median redshift of z ~ 2.4. 
We use our data to measure the Ha properties of these systems and to gauge the prevalence of active galactic 
nuclei (AGN) in these galaxies through their [Nll]/Ha ratios and Ha line widths. Removing obvious AGN, we 
find that the predicted Ha star formation rates in this diverse population are suppressed (by a factor of ~ 10) 
compared to those derived from their far-infrared luminosities. Using the AGN indicators provided by our near- 
infrared spectra, we estimate that AGN are present in at least 40% of the galaxies in our sample. To further 
investigate this, we construct a composite rest-frame spectrum for both the entire sample and for those galaxies 
which individually show no signs of nuclear activity. We find [Nll]/Ha ratios for both composite spectra which 
suggest that the energy output of the galaxies is star-formation- rather than AGN-dominated. However, we also 
find that the Ha line in the composite non-AGN spectrum is best fit with an underlying broad-line component 
with a narrow/broad flux ratio of 0.45 ± 0.20. The median Ha line width for our sample (removing obvious AGN) 
is 400 ± 70 km s" 1 (FWHM), and the typical spatial extent of the Ha emission in our narrow -band observations is 
< 4-8 kpc, which indicates a dynamical mass of 1-2 x 10 1 'M a with corresponding dynamical times of 10-20 Myr. 
Using both high-resolution imaging and spectroscopically identified velocity offsets, we find that seven of the 
far-infrared luminous galaxies have companions, suggesting that they are undergoing interactions/mergers and 
from their relative velocities we can determine a dynamical mass of 1.5 ±0.9 x 10 M Q . These measurements are 
comparable to millimeter CO estimates for the dynamical masses of these systems on similar scales, and larger 
than recent estimates of the dynamical masses of UV-selected galaxies at similar redshifts derived in an identical 
manner. Using the [Nll]/Ha index to predict abundances, we investigate the Luminosity-Metallicity relation for 
these galaxies and find that many have metallicities consistent with UV-selected high-redshift galaxies and slightly 
lower than local luminous infrared and elliptical galaxies (although we caution that our metallicity estimates have 
possible systematic uncertainties). We also compared our Ha and far-infrared luminosities with deep Chandra 
observations of a subset of our survey fields and use these data to further assess their AGN content. We conclude 
that these high-redshift, far-infrared luminous galaxies represent a population of massive, metal-rich, merging 
systems with high instantaneous star formation rates, strong dust obscuration and actively-fueled AGN which are 
likely to be the progenitors of massive local elliptical galaxies. 

Subject headings: galaxies: active — galaxies: evolution — galaxies: high-redshift — galaxies: starburst — 
submillileter 



1. INTRODUCTION 

Recent surveys in the submillimeter (sub-mm), millimeter 
(mm) and radio wave bands suggest that the star formation den- 
sity detectable by these dust-independent tracers has evolved 
strongly with redshift. Indeed, this evolution appears to out- 
strip that found using tracers which are more sensitive to dust 
obscuration, suggesting that an increasing proportion of activ- 
ity in more distant galaxies may be highly obscured (Blain et al. 
1999b). The populations resolved in these wavebands appear to 
be responsible for much of the energy density in the extragalac- 
tic far-infrared/sub-mm background (Smail et al. 2002; Cowie 
et al. 2002; Chapman et al. 2004a). However, the extreme faint- 
ness of optical counterparts to these obscured galaxies has made 
it very difficult to obtain accurate redshifts and measure intrin- 
sic properties (e.g., Simpson et al. 2004). 

The best-studied examples of the high-redshift, far-infrared 



luminous galaxy population are those identified in the sub-mm 
wave band using the SCUBA camera on the JCMT The me- 
dian redshift for this population is < z >~ 2.4 (Chapman et al. 
2003a, 2003b, 2004a), and their sub-mm fluxes suggest that 
they have bolometric luminosities > 10 12 L Q - implying that 
they are Ultraluminous Infrared Galaxies (ULIRGs). The na- 
ture of this population and their relevance to models of galaxy 
formation models and evolution are particularly important (e.g., 
Genzel et al. 2003; Baugh et al. 2004). If they are powered 
purely by star-formation then these galaxies form about half of 
the stars seen locally (Lilly et al. 1999). However, both AGN 
activity and star-formation could contribute to their immense 
far-infrared luminosities, and without further information it is 
impossible to disentangle the precise energy source (Alexander 
etal. 2003a). 

Rest-frame optical emission line properties provide a pow- 
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erful tool to investigate star formation rates (SFRs), power 
sources, and metallicity of galaxies. In particular, the Hydrogen 
Balmer emission line series is one of the primary diagnostics of 
the SFR in nearby galaxies, with the strength of Ha and its 
relative insensitivity to extinction making it the line of choice. 
The nebular recombination lines are a direct probe of the young, 
massive stellar population, since only stars with masses >10M Q 
and lifetimes < 20 Myr contribute significantly to the integrated 
ionising flux. Hence, the strength of these emission lines pro- 
vides a nearly instantaneous measure of the SFR, independent 
of the previous star formation history. Moreover, by combining 
the SFR inferred from this diagnostic line with the far-infrared 
emission (which comes from reprocessed radiation which has 
been absorbed and re-emitted by dust in the far-infrared at 
wavelengths of 10-300/im), we can gauge the prevalence of 
dust obscuration. The width of the Ha line and its intensity 
relative to other rest-frame optical lines can also give impor- 
tant information about the presence and luminosity of an AGN 
within a galaxy (Veilleux et al. 1987, 1995). 

Ha is visible in the near-infrared wave band out to z ~ 2.6 
and projects exploiting the Ha emission at these high redshifts 
have provided unique insights into the star formation properties 
of distant galaxies (e.g., Yan et al. 1999; Erb et al. 2003; Shap- 
ley et al. 2004; van Dokkum et al. 2004). With precise red- 
shifts for the far-infrared luminous population from the work 
of Chapman et al. (2004a, 2004b), we can efficiently target the 
Ha emission from the galaxies to understand the formation and 
evolution of this population and can identify the power sources, 
star formation rates, metallicities (which are accessible through 
the [Nll]/Ha [N2] index; Pettini & Pagel 2004), and masses, 
as well as more general issues such as their relation to other 
classes of high-redshift sources such as Lyman Break Galaxies 
(LBGs; Pettini et al. 2001; Erb et al. 2003; Shapley et al. 2004). 

As well as measuring the star formation rate from Ha, it is 
also possible to measure dynamics of these systems from the 
same spectra. Far-infrared luminous galaxies at z ~ 2 appear 
morphologically complex (Chapman et al. 2003b; Smail et al. 
2004). By measuring the internal dynamics of the galaxies or 
the velocity offsets between companions, we can also place lim- 
its on their masses and so test if these galaxies are truly massive 
systems (Genzel at al. 2003). 

In this paper we present the results from a near-infrared study 
of a sample of far-infrared- detected galaxies at z= 1.4-2. 7. We 
use the Ha emission line to derive their star formation rates 
and dynamics. We investigate the metallicity of these galaxies 
through their [Nll]/Ha emission and dynamics as traced by ve- 
locity structures visible in Ha emission in a subset of galaxies. 
Using Chandra data, we also compare our Ha and far-infrared 
luminosities with X-ray luminosities in several of our survey 
fields. 

In §2 we present the data reduction and results from the 
spectroscopic survey and narrow-band imaging. In §3 we dis- 
cuss the Ha properties of the far-infrared luminous galaxies 
and present a discussion of the star formation rates estimated 
from their Ha and far-infrared emission and of their dynamics, 
metallicities, and X-ray counterparts. We present our conclu- 
sions in §4. We use a cosmology with H () = 72kms _1 , flo = 0.3 
and Ao = 0.7 in which 1" corresponds to 8.2 kpc at z = 2.4. 

2. OBSERVATIONS AND ANALYSIS 



Our target sample comes from two catalogs of far-infrared 
luminous galaxies by Chapman et al. (2003a, 2004a). The 
majority of our sample comprises sub-mm detected, radio- 
identified galaxies (SMGs). These galaxies have precise po- 
sitions from the microjansky (/iJy) radio emission and are con- 
firmed to be far-infrared luminous from their detection in the 
sub-mm/mm wavebands with SCUBA or MAMBO. The radio- 
detected subset of the SMG population represents ^60%-70% 
of all SMGs brighter than > 5 mJy (Ivison et al. 2002; Chapman 
et al. 2001; Wang et al. 2004). In addition, we have also in- 
cluded a small number of optical and sub-mm faint radio galax- 
ies (OFRGs) at similar redshifts which have been proposed to 
be similarly luminous far-infrared galaxies, but with somewhat 
hotter dust temperatures resulting in them having comparatively 
faint submm fluxes (Chapman et al. 2004b; Blain et al. 2004). 
The median redshift of the combined sample is <z>= 2.4 ±0.2. 
For the purposes of this study we chose targets whose redshifts 
place Ha in spectral regions which are relatively free from 
strong atmospheric absorption and emission. 

We have explored two routes to investigate the Ha emis- 
sion from SMGs/OFRGs. First, we used narrow-band imag- 
ing to assess the Ha fluxes of galaxies in our sample - from 
both a tunable filter (NSFCAM on IRTF) and more traditional 
narrow-band filters (with the UFTI near-infrared imager on 
UKIRT) for galaxies whose redshifts serendipitously place Ha 
in the filter bandwidth. These narrow-band observations pro- 
vide the opportunity to determine the total Ha emission from 
the SMGs/OFRGs and search for any extended emission or spa- 
tial companions. 

Second, we used classical long-slit near-infrared spec- 
troscopy with the NIRSPEC spectrograph on Keck and the 
ISAAC spectrograph on the VLT. These observations allow us 
to measure precise systemic redshifts for these SMGs/OFRGs, 
with much higher reliability than that available from their rest- 
frame UV emission - which frequently show velocity shifts of 
100's of kms -1 . These precise redshifts are necessary for in- 
terferometric CO follow-up of these galaxies (Neri et al. 2003; 
Greve et al. 2004). The long-slit spectroscopic observations 
also allow us to measure the Ha luminosities and SFR(Ha) and 
to gauge the prevalence of AGN in these galaxies - through the 
detection of broad lines and extreme [Nll]/Ha flux ratios. 

These two approaches therefore provide complimentary in- 
formation on the Ha emission properties of the SMG/OFRG 
population. 

2.1. Narrow-band Imaging 
2.1.1. IRTF 

Narrow-band imaging of five targets was carried out using 
the NASA Infra-Red Telescope Facility 5 (IRTF) 3-m Telescope 
between 2003 April 28 and May 02. The observations were 
made in generally photometric conditions and ~ 1" seeing. We 
used the NSFCAM camera (Shure et al. 1993) which employs 
a 256 x 256 InSb detector at 0.15" pixel -1 to give a 38" field of 
view (which probes roughly 300 kpc at z ~ 2.4). The contin- 
uously variable tunable narrow-band filter (CVF) in NSFCAM 
provides an R = 90 passband which was tuned to the galaxy 
redshifts measured from the UV spectra from Chapman et al. 
(2003a, 2004). Shorter, matched broad-band imaging was in- 
terspersed between the narrow-band exposures to provide con- 
tinuum subtraction. 



5 The Infrared Telescope Facility is operated by the University of Hawaii under Cooperative Agreement no. NCC 5-538 with the National Aeronautics and Space 
Administration, Office of Space Science, Planetary Astronomy Program. 
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~ 500 kms" 1 between the rest-frame UV and optical emission 
lines can redshift the Ha out of the narrow-band filter. Shifts 
of the amplitude are relatively rare, and so we assume that our 
observations have sampled the bulk of the Ha emission in these 
galaxies. 

2.2. Keck Spectroscopy 

We obtained near-infrared spectra of 24 SMGs/OFRGs 
using NIRSPEC (McLean et al. 1998) on the 10-m Keck 
telescope. 7 These observations were obtained on the nights of 
2003 August 3, in non-photometric conditions; in somewhat 
better conditions on 2004 February 1; in 0.5" seeing and in 
photometric condition between 2004 April 6-9 and 2004 July 
31. In all four instances, observations were made using the 
low-resolution, longslit mode with a 4-pixel wide (0.76") slit. 
Spectra were obtained in the TT-band using the N7 or K' filters. 
In this configuration the resolution is R ~ 1500. To subtract 
sky emission, the observations were made in a standard ABBA 
sequence, where the object is nodded along the slit by 10-20". 
The total integration times are listed in Table 2 (the integration 
times of individual exposures were typically ~ 600s). One goal 
of these observations was to provide precise systemic redshifts 
of the SMGs/OFRGs to compare to those measured in the UV 
and also to provide a precise redshift for our on-going CO map- 
ping program on the Plateau de Bure interferometer (Neri et al. 
2003; Greve et al. 2004). To maximise the sample size, the ex- 
posure times were therefore kept short and so the Ha lines are 
rarely detected with sufficient signal-to-noise ratio to measure 
more than their most basic characteristics. We note that the suc- 
cess rate for detecting Ha emission in the SMGs/OFRGs in the 
near-infrared is ~ 70%. This estimate is based on the results 
from the third observing run, which was largely photometric 
since the first two observing runs were plagued by clouds. In 
cases in which there was no detection of either continuum or 
Ha emission, we have assumed that the target was not on the 
slit (as a result of an offset error or loss of guiding during the 
offset from a nearby bright star). In two cases we detected con- 
tinuum from the target but could not detect any emission lines. 
In these cases we have quoted upper limits for the Ha luminos- 
ity and SFR. 

The spectroscopic observations were reduced using the 
WMKONSPEC package in IRAF. We remap the two dimensional 
spectra using linear interpolation to rectify the spatial and spec- 
tral dimensions. After subtracting pairs of nod-positions, resid- 
ual sky features were removed in IDL using sky regions on ei- 
ther side of the object spectrum. The wavelength calibration 
from the 2003 run used the night sky lines, while for the 2004 
observing runs we used an argon arc lamp. The output pixel 
scale is 4. 3 A pixel" 1 , and the instrumental profile has a FWHM 
of 15A (measured from the widths of the sky-lines). In all line 
widths quoted in the following sections we have deconvolved 
the instrumental profile from the FWHM of the galaxies. 

The redshifts of our targets have been relatively well con- 
strained from rest-frame UV spectroscopy carried out with 
LRIS on the Keck telescope. However, again we stress that 
the rest-frame UV and optical emission lines can show velocity 
offsets of a few 1 00's km s" 1 (Pettini et al. 200 1 ; Erb et al. 2003 ; 
All of our narrow-band observations were taken based on the Shapley et al. 2003; Smail et al. in prep). With knowledge of the 
rest-frame UV redshifts and assuming that the rest-frame opti- UV redshift we were able to identify the Ha and [Nil] emission 
cal redshifts were in agreement. However, a velocity offset of in each spectrum and ensure the correct line identification. 

6 The United Kingdom Infrared Telescope is operated by the Joint Astronomy Centre on behalf of the U.K. Particle Physics and Astronomy Research Council. 

'Obtained at the W.M. Keck Observatory, which is operated as a scientific partnership among the California Institute of Technology, the University of California 
and the National Aeronautics and Space Administration. The Observatory was made possible by the generous financial support of the W.M. Keck Foundation. 



The observations were taken in a standard nine-point dither 
pattern and reduced using a running flat field of the six near- 
est temporally adjacent frames, masking bright objects before 
creating the flat-field frame. The final image was made by av- 
eraging the flat-fielded frames with a 3er clip to reject cosmic 
rays. To calibrate our data, we observed UKIRT faint photo- 
metric standards (Hawarden et al. 2001). These standards were 
observed at similar air masses and using the same instrumental 
configuration as the target galaxies. As well as providing good 
flux calibrations, these observations also allow us to calculate 
the relative throughput of the narrow-band and broad-band fil- 
ters. This allows precise subtraction of the continuum contribu- 
tion from the Ha emission in the narrow-band filter. The NS- 
FCAM CVF filter is wide enough that it contains both the Ha 
and [Nil] emission lines (with the [Nil] emission line included 
at 95% of peak transmission assuming that Ha falls at the peak 
of the filter trace). The median [Nll]/Ha ratio in our spectro- 
scopic sample is 0.25, and we therefore apply a 24% correction 
to the narrow-band flux to account for the [Nil] emission line. 
Exposure times, central wavelengths, and measured Ha fluxes 
from these observations are given in Tables 1 & 2. We illustrate 
the Ha morphology for one of the more unusual galaxies in our 
sample in Fig. 4 and discuss this in more detail in §3.1 

2.1.2. UFTI 

We obtained classical narrow-band imaging of two SMGs, 
SMMJ105226.61 and SMM J131232.31, whose redshifts 
should place Ha emission in the Br7 and H2 S 1 filters, re- 
spectively. These data were taken with the UFTI near-infrared 
camera (Roche et al. 2003) on UKIRT 6 on 2003 February 25- 
27 in < 0.7" seeing and photometric conditions. The UFTI 
camera has a 1024 x 1024 HgCdTe array and a plate scale 
of 0.091" pixel -1 , giving a field of view of 92"x92" (which 
probes roughly 0.75 Mpc in our adopted cosmology). To pro- 
vide continuum subtraction, we interspersed the narrow-band 
observations with /T-band imaging. The data were reduced 
using the ORAC-DR pipeline. The relative throughputs of the 
broadband and narrowband filters and flux calibration were de- 
termined using Two Micron All Sky Survey (2MASS) /T-band 
photometry on bright (K ~ 12) stars within the 92" field of 
view. 

Table 1: Summary of Narrow Band Imaging 
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SMM J105230.73+572209.5 


2.3692 


2.5 


12.6 


SMMJ123635.59+621424.1 


1.9721 


1.4 


9.0 


SMMJ131215.27+423900.93 


2.3396 


1.4 


7.2 


SMMJ140104.96+025223.5 


2.3396 


0.7 


9.0 


SMM J163631 .47+405556.9 


2.1544 


4.5 


14.4 


UKIRT Narrow-band 








SMM J105226.61+5721 13.0 


2.248 


5.6 


22.7 


SMMJ131232.31+423949.5 


2.166 


4.0 


15.0 



Notes: X cm denotes the central position of the narrow-band filter for 
the observations. 
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As the first two and final observing runs were non- 
photometric, we calibrate the integrated Ha fluxes onto an ab- 
solute flux scale by careful comparison of the spectra with the 
TT-band photometry of the galaxies (Smail et al. 2004). We con- 
volve the galaxy spectra with the normalised TT-band filter re- 
sponse curve to obtain the mean flux conversion factor for each 
galaxy spectrum. We also note that the atmospheric transmis- 
sion at the wavelength of redshifted Ha is ~ 1 for all targets ex- 
cept SMMJ123621.27, where the transmission is 0.75 ±0.02. 
For the third observing run we used a standard star (FS 27) to 
flux-calibrate the data. We also use the broadband photometry 
to flux-calibrate these data in order to verify that the process ap- 
plied to calibrate the previous observations is reliable and ob- 
tain a mean ratio of calibration factors of 1.05 ±0.10. This 
gives us confidence that using the photometry to flux-calibrate 
gives consistent results. All of the objects included here have 
detectable continuum emission. Where the continuum or emis- 
sion lines are faint, we binned the reduced two dimensional 
frames by a factor of two in both the spatial and spectral di- 
mension to improve the contrast prior to extraction. We show 
the extracted spectra for the whole sample in Fig. 1. 

2.3. ISAA C VLT Spectroscopy 

During the night of 2003 October 21, observations of 
SMM J221733.79 were obtained in queue mode with the 
ISAAC spectrograph on the VLT 8 (Moorwood 1997). The data 
were taken using the medium-resolution (R = 3000) grating 
and a 1" slit (with an instrumental FWHM of ~6. 7 A) and 
in 0.8" seeing. In this configuration the output pixel scale is 
1.2A pixel -1 . The total integration time was 3.0ks with the data 
taken in 4 x 750s exposures in the standard ABBA sequence 
where the object is nodded along the slit by 10-20" to achieve 
sky subtraction. We reduced the data using the ISAAC data- 
reduction pipeline, which rectifies, flat-fields, and wavelength 
calibrates the two dimensional frames. Prior to extraction, the 
data were binned by a factor of two in the spectral direction 
to boost the contrast of the object. We observed the standard 
star HIP 032007 for flux calibration. The derived Ha fluxes 
and redshifts are listed in Table 2, and the spectra are shown in 
Fig. 1. 

2.4. Spectral Analysis 

To accurately determine the redshift for each galaxy, we 
fit both the continuum level and emission lines: HaA6562.8 
and [Nll]AA6548. 1,6583.0 simultaneously with a flat contin- 
uum plus Gaussian profiles using a x 2 fit and taking into ac- 
count the greater noise in regions of strong sky emission. As 
well as fitting the narrow-line Ha and [Nil] emission, we also 
attempt to identify any underlying broad-line Ha (AGN) com- 
ponent by fitting a broad Gaussian profile at the same redshift as 
the narrow-line Ha emission but only accepting the result if the 
X 2 fit is significantly better (with >90% confidence limit) than 
with no broad-line component. This allows us to deconvolve 
the narrow-line Ha flux, which may arise from star-formation, 
from the broad-line Ha flux from an AGN. For the observa- 
tions which were flux-calibrated using the TT-band photometry 
we estimate the errors in the Ha flux and continuum levels us- 
ing both the uncertainty in the /T-band magnitude and the errors 
on the best fit to the emission line, evaluated by varying the fit 
by Ax 2 = 1. The error bars on the emission-line fluxes and 
equivalent widths are therefore conservative errors which take 

8 Based on observations collected with the ESO VLT-UT1 Antu Telescope (072 



into account the error in using the TT-band photometry and the 
uncertainty which arises from the signal-to-noise ratio in the 
data. This flux error is propagated through to the star forma- 
tion rate. We list the parameters for the fits for narrow (and 
where relevant broad) components in Table 2, along with their 
uncertainties. We note that the [Sll]AA6716,6731 doublet is 
detected in at least three individual spectra: SMM J123600.00, 
SMMJ163650.43, and SMM J163639.01. 

The ratio of [Nil] to Ha line fluxes from our spectra, as 
well as the presence of a broad component, can be used to 
identify luminous AGN in these galaxies (e.g., Veilleux et al. 
1987; Armus et al. 1989). We measure all of these observ- 
ables from our spectra and use them to flag galaxies whose 
Ha and far-infrared luminosities could be affected by emis- 
sion from an AGN, rather than being star-formation dominated 
(Table 2). Using the [Nll]/Ha emission line ratio and Ha 
line width as diagnostics, we separate our sample into three 
classes: AGN dominated, intermediate, and systems whose 
properties are consistent with star-formation (but may still 
contain some AGN). In a later section we test the reliabil- 
ity of these classes through the use of deep X-ray observa- 
tions. We classify galaxies whose near-infrared spectra show 
[Nll]/Ha<0.7 and FWHM rest <500kms- 1 as star-bursts (SB). 
Galaxies with FWHM rest between 500 and 1000km s" 1 are clas- 
sified as intermediate, while galaxies with [Nll]/Ha>0.7 and/or 
FWHM lest > 1000km s" 1 are classified as AGN. Using this clas- 
sification, ^40% (9/24) of the objects in our detected spectro- 
scopic sample have some indication of an AGN (most clas- 
sifications are based on Ha line widths, with only one AGN 
changing classification if we remove the limits on [Nll]/Ha). 
This is similar to the rate estimated from UV spectroscopy and 
high-resolution radio imaging (Chapman et al. 2003a, 2003b, 
2004c). By comparing the rest-frame UV and optical spec- 
tral classifications (Table 2), we find that the fifteen galaxies 
identified as SB in the rest-frame UV are classified as SB [9], 
AGN[3] and intermediate[3] from their rest-frame optical spec- 
tra. Likewise, the three AGN identified in the rest-frame UV 
are all classified as AGN from their rest-frame optical spectra. 
We conclude that these classifications are in reasonable agree- 
ment. However, we note that their AGN could be easily hidden 
from view in both wavebands (see e.g., §3.5). Looking at the 
far-infrared luminosities of the three spectral classes we find 
median values of 5.2 ±3.0 x 10 12 L Q for the SMGs/OFRGs 
classed as AGN, 6.6 ±0.7 x 10 12 L Q and 5.4 ± 1.4 x 1O 12 L 
for the starbursts and intermediate respectively. We thus find 
no strong evidence for a strong luminosity dependence of the 
different spectral classes. 

To determine total Ha fluxes, we have also corrected for slit 
losses based on the average /f-band light distribution of the 
galaxies. As we know the position, width, and orientation of 
the spectroscopic aperture for each galaxy, we use the TT-band 
images to calculate the fraction of the total /f-band light that 
enters the spectroscopic aperture. The galaxies are frequently 
extended and disturbed (e.g., Smail et al. 2004), and while this 
correction is uncertain, we view it as more reliable to apply this 
factor before comparing Ha to far-infrared star formation rates, 
rather than to ignore it. However, we note that using a sin- 
gle factor for correcting slit losses may contribute some scatter 
when comparing the Ha and far-infrared star formation rates. 
By careful comparison of a simulated slit for each galaxy in 
our sample with the /T-band imaging, we estimate the fraction 
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FIG. 1 . — Hq spectra for the SMGs/OFRGs in our sample. The lower spectrum in each panel is the raw spectrum (offset in flux scale for clarity). The upper spectrum is smoothed to 
the instrumental resolution. The dashed vertical lines show the positions of the HaA6562.8, [Nil] A A6548. 1,6583.0, and [Sit] A A67 17,6731, emission lines. The lower axis scale shows the 
observed wavelength, while the upper scale displays the rest-frame wavelength. The shaded areas are regions of strong sky emission. The spectra are ordered in right ascension range as 
in Table 2, and we label each galaxy with its spectral classification from §2.4. Note that we identify two Ha emission lines from SMMJ123707.21 rather than Ha + [Nil] (see §3.1). The 
final spectra are higher resolution ISAAC/VLT spectra (Table 2). 



of the flux entering the slit compared to the total K-band flux 
of the galaxy, /, and obtain a mean value of / = 0.62 ± 0.06. 
We assume that the equivalent width of the Ha emission line 
is constant across the galaxy and so only those galaxies with 
spectra taken in photometric conditions, whose flux calibration 



was carried out using the standard star observations, require this 
correction to their Ha flux. 

Finally, we note that there are two repeated observations of 
targets between the spectroscopic and narrow-band imaging 
samples: SMM J105230.60 and SMM J123635.59 — as well 
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as narrow-band imaging of SMM J140104.96 for which an Ha 
spectrum exists in the literature (0.56 ±0.08 x 10" 18 Wm" 2 ; Ivi- 
son et al. 2000). The agreement between the Ha fluxes from 
the two techniques is very good (Table 2), with a median ratio 
of spectroscopic to imaging fluxes of 0.96 ±0.04. On the ba- 
sis of this good agreement we feel confident in combining the 
spectroscopic and narrow-band imaging data to discuss the Ha 
luminosities of SMGs/OFRGs in the following sections. 



2.5. Star Formation Rates 

For solar abundances and adopting a Salpeter initial mass 
function (IMF), the conversion between Ha flux and SFR is 
SFR(M Q yr 1 ) = 7.9 x 10" 35 L(Ha) W (Kennicutt et al. 1998). 
This calibration assumes that all of the ionising photons are 
reprocessed into nebular lines (i.e., they are neither absorbed 
by dust before they can ionise the gas, nor do they escape the 
galaxy). 
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Table 2: Summary of Results 
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2.4662[5] 


3.4 ±0.6 


434 ± 25 


137 ±40 


125 ±20 


< 0.05 


SMMJ123635.59+621424.1 


1.8 


2.0150[2] 


4.2 ±0.3 


240 ±33 


45 ±10 


130 ±30 


0.67 ±0.27 








6.9 ±0.9 


1623 ±213 


73 ± 10 






OJVLM J IzJ /U/.Zl+OZl4Uo. 1 


J.D 




1 £ 4_ o < 
i .o ± u.j 


_i4o zt 4U 


9< 4- Q 


CQ 4_ 9/1 

ftft ± Z4 




SMM T 1237 11 32+621331 3 


2.4 


1 9958T41 


0.4 ±0.3 


112± 18 


10 ±5 


16±9 


< 0.05 


SMMJ123711.98+621325.7 


2.4 




< 0.25 




< 5 


< 6 




SMM T131905 60+493946 


2.4 




< 0.3 




< 6 


- o 




dht 1 a i 90S 1,Aj-A1A\ AA A 


9 zl 


9 9379N SI 
Z.Zj /Z[ 1 oj 


U.V H: U.Z 


AA9. 4- fi(\ 


94 4- ^ 
Z4 it J 


99 4_ i o 
J-X 1ft 


1 90 4- 60 
1 .ZU It U.OU 








0.9 ±0.8 


959 ± 100 


26 ±5 






SMM J131208. 82+424129.1 
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0.22 ±0.08 
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2.4 


1.4880[6] 
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0.30 ±0.20 
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4.8 
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1236 ±200 






SMMJ163655.80+405914.0 


2.4 
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2.610[5] 
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15.0 
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Note the value given in the [] z column is the error in the last decimal place. The Ha flux given in column 4 is the narrow line 
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Ha fluxes are the observed flux (not corrected for slit flux losses or atmospheric transmission), but the SFR's have had both corrections applied (where applicable). 
The Ha star-formation rates are uncorrected for extinction. Ha classifications are described in §2.4 (SB=star burst; int=intermediate, AGN=AGN). 
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Rest-frame Optical Spectra of SCUBA Galaxies 




Arc Seconds wavelength (fim) 



Fig. 2.— Left: Ccombined B, V and /-band H ST ACS observations of SMM J123635.59+621424.1 with the slit position overlaid. Right: The 
position- velocity diagram around the Ha line from the NIRSPEC near-infrared spectrum of the galaxy. This shows an apparent ~100kms~ 
velocity gradient in the Ha emission line and a somewhat larger velocity gradient, ~150kms~', in [Nil]. The central component of the galaxy 
(labeled Al) has an [Nll]/Ha ratio of 0.60 ±0.15 and hosts a bright, unresolved radio source (Chapman et al. 2004a); both features suggest that 
it contains an AGN. Al is offset from A2 by ~150kms~' in velocity and around 3kpc in projection on the sky. A2 has an [Nll]/Ha ratio of 
0.2 ±0.1, suggesting that this is a star-forming knot or close companion. Here 1" corresponds to 8.4 kpc at the redshift of this galaxy and the 
image has been rotated to match the position angle of the spectrum. 



We also have a second SFR indicator for our sample of 
SMGs/OFRGs - their far-infrared (FIR) luminosities. A sig- 
nificant fraction of the bolometric luminosity of the most active 
dusty galaxies is absorbed by interstellar dust and re-emitted 
in the thermal IR, at wavelengths 10-300/im. If young stars 
dominate the radiation output in the UV-visible wave band and 
the dust opacity is high everywhere, then the far-infrared lumi- 
nosity measures the bolometric luminosity of the star-burst and 
this in turn provides an excellent tracer of the SFR of the galaxy. 
Adopting the models of Leitherer et al. (1995) for continuous 
bursts of age 10-100 Myr and using the same IMF as in the Ha 
calculation yields SFR(M Q yr" 1 ) = 4.5 x 10" 37 L(FIR)W (Ken- 
nicutt et al. 1998). We stress that this relation only holds if the 
age of the star-burst is less than 100 Myr. 

We list the SFR and Ha flux measurements from our sample 
in Table 2. We also give the far-infrared luminosities of these 
galaxies from Chapman et al. (2003a, 2004b). These values are 
derived from fitting model spectral energy distributions (SEDs) 
to the observed 850 /im and 1 .4-GHz fluxes of the galaxies at 
their known redshifts, assuming that the local far-infrared-radio 
correlation holds (Condon et al. 1991; Garrett 2002). 

2.6. Notes on Individual Galaxies 

Observations of six of the galaxies in our sample are partic- 
ularly noteworthy and so we discuss these in more detail here. 

SMMJ123635.59+621424.1 — Dawson et al. (2003) identi- 
fied this object as a z=2.015 spiral galaxy, but their near- 
infrared spectroscopy around the Ha emission line shows 
a 2500±250 kms" 1 broad-line Ha component and [Nll]/Ha 
emission line flux ratio 0.45 ±0.1. The presence of the broad- 
line component and detection of hard (V = 0.3) X-ray emis- 
sion from Chandra imaging indicate an obscured type II AGN 
(Dawson et al. 2003). Our spectroscopic and narrow-band ob- 
servations of this galaxy produce comparable Ha fluxes, which, 
along with the compact morphology of the galaxy in the IRTF 
narrow-band image, suggests that the spectroscopic slit is sam- 
pling the bulk of the emission from this system (the slightly 
lower flux in the narrow-band observation arises because the 



narrow-filter only samples the narrow-line Ha; <~ 20% of the 
broad-line flux is missed). Our spectroscopic observations of 
this target show clear velocity structure in the Ha and [Nil] 
emission lines (Fig. 2). We aligned the slit with the center of 
the galaxy and the bright knot (labeled A2 in Fig. 2), and our 
resulting Keck spectroscopic observations indicate either a ve- 
locity shear or rotation in the Ha emission line with an ampli- 
tude of <~ 100 km s" 1 within ~ 2" (17 kpc in projection) along 
the slit (Fig. 2). By collapsing down the Ha and [Nil] emis- 
sion lines in the spectral direction, we find an offset of 0.3" 
(3 kpc) along the slit between the maximum intensities of these 
two lines, as well as an associated variation in the [Nll]/Ha ra- 
tio from ~ 0.2 ±0.1 up to <~ 0.60 ±0.15 (which is consistent 
with the results of Dawson et al. 2003). We also find a similar 
broad-line component to the Ha emission, confirming the pres- 
ence of an AGN. The high-resolution Hubble Space Telescope 9 
(HST) ACS image of this galaxy (Fig. 2), from the GOODS 
imaging of this region (Giavalisco et al. 2004; Dawson et al. 
2003; Smail et al. 2004), shows an apparently face-on spiral 
galaxy with a bright nucleus (Al) and a prominent companion 
or knot in one of the spiral arms (A2); the optical extent of the 
galaxy is 13 kpc (1.5"). The separation between the nucleus 
(Al) and knot (A2) in the image is <~ 0.5" (5 kpc) - compara- 
ble to the apparent offset in the emission-line peaks. The high- 
resolution 1 .4-GHz MERLIN map of this galaxy in Chapman 
et al. (2004c) shows an unresolved radio source coincident with 
the nuclear component, which is also where the [Nll]/Ha is the 
strongest (<~ 0.6), providing further support for the classifica- 
tion of this component as an AGN. The face-on aspect of this 
system, combined with the modest velocity difference between 
Al and A2, suggests to us that the latter may be a dynamically- 
distinct component (rather than a star-forming knot inside a spi- 
ral arm), an interaction which has prompted the activity we see. 
We note that it is possible for intensity gradients between sep- 
arate components to mimic velocity gradients as a result of the 
way in which long-slit spectroscopy mixes spatial and spectral 
domains. However, the spatial offset between Al and A2 in the 
dispersion direction is ~ 0.4", corresponding to 10A, which is 
much less than the apparent velocity gradient. We therefore 



Swinbank et al. 



9 




Arc Seconds 
Center: R.A. 14 01 04.91 Dec +02 52 24.2 



Fig. 3. — Three views of SMM J140104.9: broadband K image from IRTF (left); Ha continuum-subtracted, narrow-band image from IRTF 
(middle) and the /f5r^702-band of the system which illustrates its very complex morphology (right). We identify in the HST image the three main 
components of this galaxy as defined by Ivison et al. (2001), and in each panel we also mark the position of the peak of component Jl. It is clear 
that the optical-near-infrared colors and Ha surface brightness of component J2 are much different from those of Jin and Jlse - ruling out the 
suggestion that J2 represents a lensed counter-image of Jin or Jlse (Downes & Solomon 2003). Moreover, the Ha emission extending south from 
Jin and wrapping around Jl traces the morphology of the clumps visible in the R-bwad, Jlse, and has no significant contribution from the smooth 
component, Jl. We discuss the interpretation of this result in more detail in §3.1 



suggest that the velocity offset most likely arises from motions 
within the galaxy, rather than as an artifact of the observation. 

SMM J '123707. 21 +621408.1 — We observed this target with a 
position angle such that the NIRSPEC slit passes through the 
two components shown in the HST ACS imaging in Smail et al. 
(2004). This galaxy consists of a red component and a much 
bluer object separated by 0.2" (1.7kpc; labeled Bl and B2 in 
Fig 5). The resulting NIRSPEC spectrum shows two strong, 
spatially-extended lines separated by ~ 600 km s" 1 . We iden- 
tify both of these emission lines as Ha (rather than Ha and 
[Nil]) from two separate components for two reasons: (1) the 
spatial separation between the two components in the spectra 
is matched almost exactly by the separation in the HST ACS 
image (Fig. 5); and (2) if the higher velocity component was 
identified as strong [Nll]A6583 then the line ratio would be 
[Nll]/Ha ~ 1, and yet we see no signs of the [Sll] emission 
line, which has a ratio of [Sll]/[Nll]~ 0.5 for AGN (Ferland & 
Ostenbrock 1986). 



ratio shows marginal evidence for variation from 0.45 ± 0. 10 to 
0.35 ±0.10 between CI and C2 and also mimics the velocity 
offset between CI and C2. These modest [Nll]/Ha emission 
line ratios indicate that both components are likely to be star- 
burst, rather than AGN-, dominated. 

d(kpc) 

-8.3 8.3 



2.136F 




SMMJ123711.32+621331.3 — This target was known to con- 
sist of two radio sources, both of which may contribute to 
the far-infrared emission. Aligning the slit along both com- 
ponents, we detect only faint continuum and no lines from 
the UV-identified component, SMMJ12371 1.98+621325.7, for 
which Chapman et al. (2004a) measured the redshift of the 
system. However, we do detect Ha emission from the sec- 
ond radio source, SMM J12371 1.32+621331.3. This allowed 
us to derive a redshift of 1.9958 ± 0.0004, giving an off- 
set of 400±50kms _1 and 8" (^70kpc) in projection from 
SMMJ123711.98+621325.7. 

RGJ 131236.05+424044.1 — The ground-based ZT-band imag- 
ing of this galaxy from Smail et al. (2004) shows a bright nu- 
cleus surrounded by a diffuse halo approximately 2" (16kpc) 
in extent. The two components seen in the near-infrared spec- 
trum are separated by ~ 185 ±45 km s" 1 in velocity and 0.4" 
(3.4kpc) in projection. The two-dimensional near-infrared 
spectrum suggests that one of the sources (labeled CI) displays 
signs of a velocity shear in Ha across <~ 0.5". The [Nll]/Ha 

'Based on observations made with the NASA/ESA Hubble Space Telescope 
Astronomy, Inc., under NASA contract NAS5-26555 



2.123 f. ^ , ,,,,,,, 7 

-1 1 

Arc Seconds 

FIG. 4. - Two Dimensional near-infrared spectrum of RGJ J131236.05, show- 
ing the velocity structure seen in Ha and [Nil]. We identify two components 
(labeled CI and C2), offset by 185 ±40kms~ 1 and 3.4 kpc in projection along 
the slit. 

SMM J1401 04.96+025223.5— This z = 2.56 sub-mm selected 
galaxy (SMM J1401 1+0252) is discussed in detail by Ivison et 
al. (2000, 2001). It lies in the field of the z = 0.25 cluster A 1835 
and is expected to be amplified by a factor of 2.75 by the fore- 
ground cluster potential. The morphology of this galaxy is com- 
plex (Fig. 3). Ivison et al. (2001) identify three main compo- 
nents: Jl, a blue relatively smooth and regular object; J2, a 
bluer and more compact object about 2" to the west of Jl; and 
Jin, an extremely red, diffuse structure which extends 2-3" to 
the north of J 1 . 

The interpretation of this multi-component system is con- 

which is operated by STSCI for the Assosciate of Universities for Research in 
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Rest-frame Optical Spectra of SCUBA Galaxies 



Ha Ha 
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Fig. 5. — Left: Position-velocity diagram around the Ha line from the near-infrared spectrum of SMM J123707.21. Right: Combined B-, V- and 
/-band HST ACS observations of this galaxy with the slit position overlaid. Middle: Intensity distribution of the Ha emission collapsed over the 
central 200km s -1 of each of the two Ha emission lines. This shows an apparent spatial offset of 0.2" (and 600 km s -1 velocity offset) between the 
two galaxies. This matches the apparent separation of the central component of the galaxy (labeled Bl) from the much redder second component 
(labeled B2) which are offset by 2 kpc in projection. (For a color version of this image, see Smail et al. 2004). 



tentious. Ivison et al. (2001) suggest that the gas reservoir and 
source of the far-infrared emission from this galaxy reside in 
Jin - with J2 and Jl simply being less-obscured components 
within the same system. However, Downes & Solomon (2003) 
suggest that the clumpy features to the south-west of Jl (called 
Jlse), as well as Jin and J2, are all images of a single, highly- 
amplified background sub-mm source - which is being lensed 
by both the regular component of Jl (which they identify as a 
foreground (z <~ 0.20) galaxy) and the cluster potential. 

Our IRTF continuum-subtracted Ha image provides a power- 
ful tool for testing these competing suggestions (see also Tecza 
et al. 2004). We compare the morphology of the Ha emis- 
sion with the Rjq2- and TT-band images in Fig. 3. It is clear 
that the bulk of the emission traces components Jl and Jlse, 
with the regular J 1 and compact J2 being undetected. This im- 
mediately rules out the possibility that J2 is a lensed counter- 
image of Jin or Jlse. However, the absence of Ha emission 
from the smooth component of Jl is a concern - suggesting 
that the claim by Downes & Solomon (2003) that it is a fore- 
ground lens may be correct. To further test this, we have re- 
turned to the optical spectrum of Jl presented by Barger et al. 
(1999) and claimed by them to represent a classical Lyman- 
break galaxy. More careful study of this spectrum leads us to 
re-examine this interpretation: there are strong absorption fea- 
tures at 4745, 4784, 4912, 4957, 5374, and 6070A, which are 
unidentifiable at the redshift claimed by Barger et al., z = 2.55, 
but correspond exactly to Balmer H£, Hrj, Ca H&K, G-band 
and H/3 absorption at z = 0.248. This is an unfortunate conjunc- 
tion - Jl is a ~ 0.1L*, post-starburst member of the A 1835 
cluster. The absence of any strong emission lines in this galaxy 
and broadband colours which are much bluer than the typical 
passive E/S0 cluster member mean that its nature is not imme- 
diately obvious from either the published spectrum or the true 
colour images of the cluster (Ivison et al. 2000, 2001). 

How does this effect the interpretations of SMM J140104.96 
by Downes & Solomon (2003) and Ivison et al. (2001)? We 
confirm that Jl is a foreground lens as stated by Downes & 
Solomon (2003). However, our observations disprove their cen- 
tral claim that Jin, Jlse, and J2 are highly-amplified multiple 
images of an intrinsically low-luminosity sub-mm source. In 
fact, Jln/Jlse and J2 are probably single images of three back- 
ground galaxies at z = 2.56 (see Tecza et al. 2004), with the 
Ha emission from this system arising entirely from very red 



Jln/Jlse - which is also the site of a massive gas reservoir and 
hence most likely the far-infrared source (Frayer et al. 2003; 
Ivison et al. 2001). The addition of Jl to the foreground lens 
model does increase the estimate of the area-averaged amplifi- 
cation for Jln/Jse from 2.5 to ~ 5; however, this does not sig- 
nificantly alter any of the conclusions in Ivison et al. (2001). 

SMMJ163650.43+405734.5 — We also observed 
SMMJ163650.43 (ELAIS-N2 850.4, Smail et al. 2003) us- 
ing NIRSPEC for a total of 4.8 ks at two orthogonal position 
angles. The detailed kinematics of this complex merging sys- 
tem using three dimensional near-infrared spectroscopy will be 
discussed in Swinbank et al. 2005. However, we note that the 
broad Ha component is apparently offset from the narrow Ha 
component by ~ 300kms _1 . The redshift quoted in Table 2 is 
that of the narrow line Ha. 



2.7. Ha Properties ofSMGs 

We show near-infrared spectra of the 23 detected galaxies 
from our Keck and VLT observations in Fig. 1 marked with the 
expected redshifts of lines based on the best-fit Ha redshift (Ta- 
ble 2). As these were short exposures, primarily meant to derive 
redshifts, the signal-to-noise ratio on the individual galaxies is 
generally modest. To overcome this, we have also combined all 
of the spectra to provide a composite near-infrared spectrum for 
a representative far-infrared luminous galaxy at z ~ 2.4. 

We create the composite spectrum by de-redshifting and 
summing all of the spectra (normalised by Ha flux) for our 
galaxies (we note that stacking the spectra based on their in- 
dividual signal-to-noise ratio or an unweighted stack does not 
alter any of the conclusions below). We also derive a com- 
posite spectrum for those galaxies which individually show no 
signs of an AGN (i.e., those with small [Nll]/Ha ratios and line 
widths). The resulting composite spectra are shown in Fig. 6. 
The rest-frame composite spectrum from the entire sample is 
best fit (> 99% confidence level) with an underlying broad- 
line region with a narrow-to-broad line flux ratio of 0.6 ±0.1 
and FWHM rest of 1300 ±210kms- 1 for the broad-line Ha and 
325±30kms _1 for the narrow-line Ha. The average [Nll]/Ha 
for the entire sample is 0.42 ±0.05 which is indicative of star 
formation. 
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Fig. 6. - Rest-frame composite spectrum of all of the galaxies in our sample 
(top) and well as the composite from the (individually spectroscopically clas- 
sified) star-forming galaxies (bottom) (we have not included SMM J123707.21 
in either composite as we believe that the spectra display two Ho components 
coming from two interacting galaxies rather than Ha and [Nil]). The insets 
show the region around the Ha emission line with the best fit to the emission 
line using four Gaussian profiles overlaid. We also show the broad line com- 
ponent of the best-fit. The average [Nll]/Ha ratio is 0.42 ± 0.05 for the entire 
sample and 0.19 ±0.05 for the star-forming galaxies. Both are consistent with 
star-forming galaxies rather than an AGN. We also detect the [SlI]AA6716,6731 
emission doublet (marginally detected in the SB sample) and the (unresolved) 
stellar NaAA5889. 95,5895. 92 absorption doublet. By fitting the Na absorp- 
tion and [SlI] doublets, we find no velocity offset from Ha. The [Sll]/Ha 
ratio can be used to classify the SMGs/ORFGs optical emission line properties 
as LINER- or Hll-region-like. While the SMGs/OFRGs which are individu- 
ally classified as star-forming galaxies show no signs of an underlying broad 
Ha line, the composite from this subsample of galaxies is best fit with a nar- 
row/broad Ha emission line ratio of 0.45 ± 0.20. 

The [Sll]/Ha ratio can be used to classify the spectral type 
of galaxies (Veilleux et al. 1995). The wavelengths of the 
[Sll]AA6716,6731 lines in the composite spectra are 6722 ±6 
and 6729 ± 5A - indicating no detectable velocity offset be- 
tween the Ha and [Sll] emission lines. From the strengths of 
the lines we estimate a ratio of [Sll]/Ha = 0.10 ±0.04 - plac- 
ing the composite SMG/OFRG within either the LINER or Hll 
region of the classification space (Veilleux et al. 1987, 1995). 
This is similar to the typical optical spectral classification of 
local ULIRGs, for which mid-infrared observations (Lutz et 
al. 1999) and spatially-resolved spectroscopy (Heckman et al. 
1990) both suggest that star-formation is the dominant power 
source (see Lutz et al. 1999). 

The composite spectrum also shows an absorption feature at 
a wavelength close to that expected for the Na doublet which 
is seen from cool stars and from the warm neutral (T~ 10 4 K) 
phase of the interstellar medium (ISM; Phillips 1993). This 
has been used to map the velocity structure of large-scale out- 
flows from nearby ULIRGs (Heckman et al. 2000; Rupke et 
al. 2002). For a galaxy whose light is dominated by a young 
star-burst, as we believe is the case for these SMGs/OFRGs, 



the warm ISM is expected to be the primary source of this ab- 
sorption feature. We fit the Na absorption line with two Gaus- 
sian profiles with a fixed separation and variable width and de- 
rive a central wavelength of 5889±5A, which is entirely con- 
sistent with the rest-frame wavelength of the resolved doublet 
(AA5889.95,5895.92A). We thus place a limit of <150kms" 1 
on the possible velocity offset of this feature from the Ha red- 
shift. The FWHM rest of the Na lines is 292 ± 192kms" 1 , which 
is also consistent with the FWHM of the Ha emission in the 
rest-frame composite. 

The composite spectrum thus shows no signs of an offset in 
velocity between Na or [Sll] and Ha. This is slightly surpris- 
ing since many local luminous/ultra-luminous infra-red galax- 
ies display velocity offsets of several 100's kms -1 between 
these lines (Heckman et al. 2000; Rupke et al. 2002). We de- 
fer a detailed discussion of the velocity offsets between the UV 
ISM and (star-burst) nebular emission lines in the individual 
galaxies to a later paper (Smail et al. in prep). 

Turning to the rest-frame composite from those galaxies 
which individually show no signs of an AGN component (i.e. 
those with low [Nll]/Ha emission line ratios and small Ha 
line widths), we fit this spectrum with both [Nll]AA6548,6583 
emission lines and also attempt to fit a broad line Ha compo- 
nent. When the broad component is included in the fit, the total 
X 2 is marginally better than when the broad component is ex- 
cluded, with a change in the total x 2 = 4 which corresponds 
to ~ 2(7 (or ~ 87% confidence). The resulting best-fit model 
for the SB composite has an underlying broad-line Ha com- 
ponent with a ratio of broad-line/narrow-line Ha = 0.45 ± 0.20 
and a broad-line FWHM rest = 890 ±2 10 kms -1 , suggesting that 
even the SMGs/OFRGs which are identified individually as 
star-bursts may contain at least some level of underlying non- 
thermal activity. Nevertheless, the [Nll]/Ha emission line ratio 
is 0.19 ±0.05, which, along with the limit of [Sll]/Ha < 0.14, 
still suggests that the energy output is star-formation-, rather 
than AGN- dominated. 

2.8. Kinematics of Ha emission 

The narrow Ha emission line of the SB composite has a 
rest-frame FWHM, est , after correcting for the instrumental res- 
olution of the observations, of FWHM lest = 400 ±70 kms" 1 
- in agreement with the average Ha line width from the 
SMGs/OFRGs which make up this composite spectrum (350 ± 
50kms _1 ). This is somewhat larger than the Ha line widths of 
UV-selected galaxies at z ~ 2 identified by Erb et al. (2003), 
who derive a mean FWHM rest = 242 ± 65 km s" 1 for their sam- 
ple. We compare the distribution of line widths for these two 
populations (as a function of their Ha luminosities) in Fig. 7. 
We see that the SMGs/OFRGs are typically 6 times brighter in 
Ha (before any correction for reddening) than the UV-selected 
population at their epoch. More interestingly, the emission line 
widths of the SMGs/OFRGs are on average ~ 50% larger than 
those measured for the UV-selected systems, although the two 
distributions overlap substantially. This difference could re- 
flect either: (1) different halo masses, (2) different dynamical 
states, (3) a contribution from a broad AGN component in some 
SMGs/OFRGs, or (4) starburst-driven superwinds contributing 
to the line widths. The composite spectra from §3.2 certainly 
suggest that an unidentified broad component to the Ha emis- 
sion may be present in some SMGs/OFRGs - however, as we 
show below, we do not believe that this is the chief cause of the 
difference in the line widths of the two populations. 

We first compare our Ha-derived line widths with those ob- 
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FIG . 7 . — Left: Comparison of the Ha luminosity vs FWHMrest for the narrow Ha components (deconvolved for instrumental resolution) in our sample compared 
to those found in UV-selected galaxies by Erb et al. (2003; we have conservatively assumed a factor of two uncertainty in the Ha luminosities for this sample). We 
also include those galaxies which show single broad lines (no narrow component). The crosses mark the positions of the median luminosity and FWHM rest for the 
two populations. The line widths in the SMGs/OFRGs have a median of 330kms and a long tail out to ~ 600 kms~' . In comparison to the UV-selected galaxies 
studied by Erb et al. (2003) (which have a median of ~ 210kms~'), SMGs/OFRGs appear to be typically more massive systems, with a wider distribution of line 
widths (although some of the broadest lines may come from unresolved companions or AGN). Right: Correlation between rest-frame optical emission line width 
and rest-frame V-band luminosity (the Faber-Jackson relation) for the SMGs/OFRGs compared to the LBGs and DRGs from Pettini et al. (2001) and van Dokkum 
et al. (2004). The crosses mark the positions of the median luminosity and <7 res t for the SMGs/OFRGs and LBGs. We show the local Faber-Jackson relation by the 
dotted line (assuming V -R=0.3) and also the same line offset by a further 2.5 magnitudes in order to pass through the median of the SMG sample (solid line). 



tained from the dynamics of cold gas as traced by interferomet- 
ric maps of the CO distributions in a small number of SMGs 
(Frayer et al. 1999; Neri et al. 2003; Genzel et al. 2003). The 
sample of five galaxies compiled by Neri et al. (2003) have a 
FWHM rest for the CO of 420±35 km s _1 and a mean ratio of Ha 
to CO FWHM consistent with unity (1.27 ±0.24) for the four 
galaxies in common with our sample. The CO line widths are 
unaffected by the presence of any AGN or large-scale winds, 
and so this suggests that the AGN or superwind contribution to 
the Ha line widths for SMGs/OFRGs may not be responsible 
for the comparatively large Ha line widths. 

We can therefore combine the emission-line widths for our 
sample with the typical physical extent of the Ha emission from 
our narrow-band imaging of the SMGs/OFRGs to place limits 
on their masses. The spatial extent of the Ha in the galaxies has 
a wide distribution: < 0.5-1.0" (corrected for seeing) or < 4- 
8 kpc (c.f. Smail et al. 2004). Assuming that the Ha emission 
arises from virialised clouds in the galaxy's potential well, we 
estimate a typical mass of 1-2 x 10 1 1 M Q for the SMGs/OFRGs 
in our sample (Erb et al. 2003), with corresponding dynamical 
times of 10-20 Myr. Using the limits on the spatial extent of the 
CO emission, Neri et al. (2003) determine a median dynamical 
mass of ~ 6 x 10 10 M Q assuming that the CO gas is in bound 
orbits. Thus, the masses derived from the dynamics of the cold 
gas in a small sample of these galaxies support those estimated 
from the emission-line kinematics. 

The masses of UV-selected galaxies at z ~ 2 derived in 
an identical manner by Erb et al. (2003) have a median of 
3.3 ± 1.1 x 10 10 M Q , around 5 times lower than our estimates 
for the SMGs/OFRGs. This is due to a combination of smaller 
estimated sizes and lower Ha line widths. Clearly both of these 
estimates have large systematic errors, yet they are suggestive 
of a real difference in the characteristic masses (or dynamical 



states) of restframe UV- and far-infrared selected galaxies at 
z>2. 

We can also obtain an independent measure of the masses of 
the SMGs/OFRGs from the seven galaxies which show multi- 
ple components and for which we have estimates of the velocity 
differences between the components. Assuming random orien- 
tations on the sky and that the components are bound/merging, 
the masses of the systems are ~ 1.5 ±0.9 x 1O U M0. This is 
comparable to the earlier estimates and gives us confidence that 
the SMGs/OFRGs are indeed massive galaxies. 

2.9. Magnitude - Line Width correlations 

Correlations between rest-frame luminosity and kinematics 
(as measured by the widths of the emission lines) at high red- 
shift have had only limited success. Pettini et al. (2001) and 
van Dokkum et al. (2004) have attempted to study the Tully- 
Fisher-like or Faber-Jackson-like (Tully & Fisher 1977; Faber 
& Jackson 1976) correlations for LBGs and distant red galaxies 
(DRGs) at z>2; however, no correlations have been found over 
the range FWHM les t ~ 120- 300kms _1 . This may be due in 
part to the small sample sizes involved. To search for such a 
correlation, the SMGs for which we have well-defined narrow- 
line Ha line width measurements can be added to this sample. 
We compute the de-reddened rest-frame V-band magnitudes for 
the SMGs in our sample by using HYPER-Z (Bolzonella et al. 
2000) to compute the best-fit SED to the observed UK photom- 
etry from Smail et al. (2004). We also convert FWHM rest to a 
by assuming FWHM= 2.35 x a. 

Fig. 7 shows the resulting luminosity-line width correlation 
for the combined sample of SMGs, LBGs, and DRGs (the latter 
two from Pettini et al. (2001) and van Dokkum et al. (2004) are 
corrected for reddening using their estimates of Ay). We also 
overlay the local Faber-Jackson relation from J0rgensen et al. 
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(1995, 1999; assuming V-R = 0.3), and fit the zero-point of 
the same correlation so that it passes through the median of the 
SMG sample - resulting in an offset of 2.5 mag. This offset is 
comparable to the Af ac [ e V found in Smail et al. (2004) and sug- 
gests that the descendants of these high redshift populations are 
likely to lie on or around the local Faber-Jackson relation. 

2.10. SFR Comparisons 

Next we compare the far-infrared and Ha luminosities of the 
galaxies in our sample to investigate the influence of AGN and 
star formation power-sources and the possible effects of dust 
extinction on the Ha emission from this population. 

Fig. 8 compares the Ha and far-infrared luminosities of the 
galaxies in our sample. We have also included measurements 
from the literature from five well-studied dusty, high-redshift 
galaxies: SMMJ02399-0136 (Ivison et al. 1998, 2000), 
SMMJ17142+5016 (Smail et al. 2003), ERO J 164502+4626 
(Dey et al. 1999), SMMJ0443 1+0210 (Frayeret al. 2003), and 
SMM J16359+6612 (Kneib et al. 2004). 

The star formation rate derived from the far-infrared and 
Ha should be correlated if the effects of dust and any con- 
tributions from AGN are uniform across the sample. The 
SMGs/OFRGs in our sample show only a weak correlation 
between SFR(FIR) and SFR(Ha), with comparable scatter in 
the Ha-derived SFR and that estimated from the far-infrared 
(A(FIR)/FIR=0.30±0.18 vs. A(Ha)/Ha = 0.37 ±0.20). Re- 
assuringly, the two galaxies for which we failed to obtain Ha 
detections in good conditions are also two of the least luminous 
galaxies when ranked on their far-infrared emission. 

We have also included in Fig. 8 samples of similar and less 
luminous far-infrared galaxies from surveys of the z < 1 uni- 
verse. These come from ISOCAM-selected galaxy surveys 
(Franceschini et al. 2003; Flores et al. 2004) and studies of lo- 
cal Very- and Ultra-Luminous infrared galaxies by Poggianti & 
Wu (2000) and Dopita et al. (2002). Compared to these sam- 
ples, we see a similar wide dispersion in the ratio of far-infrared 
to Ha luminosities across a factor of nearly 100 in far-infrared 
luminosity. This is suggestive of a similar diversity in the en- 
ergy sources and obscuration for galaxies with far-infrared lu- 
minosities from 10 10 -> 1O 12 L . In support of this, we note that 
the distribution of Ha equivalent widths (EW, Table 2) for the 
SMG sample peaks at about ~ 20-40A and has a long tail out 
to > 100A, with a median EW(Ha) = 75 ± 25 A. The shape of 
the distribution is very similar to that seen for the rest-frame 
Ha EWs of local ULIRGs from Veilleux et al. (1999), which 
have a median EW(Ha) = 73 ± 8A suggesting that the unob- 
scured/partially observed mix of emission-line gas and stellar 
continuum is comparable to local ULIRGs (although we note 
that the projected size of the spectroscopy apertures in the local 
ULIRGs is less than 1 kpc, whereas for the SMGs/OFRGs the 
projected size is ~ lOkpc; thus some caution should be taken 
when comparing the two samples). 

Our SMG/OFRG sample is intrinsically luminous in the far- 
infrared, but their Ha flux (and estimated SFR) is much less 
than expected with a median SFR(Ha) compared to SFR(IR) 
of 94 ±20 and 1380 ± 19OM0yr _1 , respectively. Overall, the 
Ha SFRs are suppressed by at least a factor of ten relative to 
that suggested by the far-infrared. In comparison to other far- 
infrared-selected samples, we see that the SMGs/OFRGs ex- 
tend the trend for proportionally less Ha luminosity in more far- 
infrared luminous galaxies. As we have shown there is no de- 
tectable difference between the far-infrared luminosities of the 
starburst- and AGN-classed galaxies in our sample - suggest- 



ing that this declining ratio of Ha to far-infrared emission is un- 
likely to be caused by an increasing AGN contribution. Hence, 
we attribute the variation to copious and increasing amounts of 
dust enshrouding the galaxy and extinguishing Ha emission. 

Unfortunately, the wavelength coverage of our spectroscopic 
observations does not extend to H/3, and thus the reddening 
in these galaxies cannot be estimated directly from the Balmer 
decrement to confirm this suggestion, although attempts at esti- 
mating the reddening in this manner are underway for a sub- 
sample of galaxies. The only target with an H/3 measure- 
ment in the literature is SMM J123707.21+621408.1 (Simpson 
et al. 2004), with an H/3 flux of 2.1 ±0.9 x 10- 20 Wnr 2 . The 
Ha/H/3 flux ratio is 8±6, which corresponds to a reddening 
of A v = 1.4 ± 1.0 (Calzetti et al. 1994); the large uncertain- 
ties arise from the modest signal-to-noise ratio detections of 
both Ha and H/3 emission lines. For the rest of the sample, 
we have attempted to derive the reddening for these galaxies 
from their broadband optical/near-infrared colours (Smail et al. 
2004). These are derived by using the HYPER-Z photometric 
redshift code (Bolzonella et al. 2000) to fit young continuous 
star formation models with variable reddening and age to the 
galaxy photometry at their known redshifts. We find a me- 
dian reddening of Ay = 3.0 ± 1.0. Accounting for the contri- 
bution of Ha to the /f-band photometry may decrease this es- 
timate slightly but the reddening estimated from the continuum 
colours indicates substantial extinction at the wavelength of Ha 
(assuming that the continuum and line emission arise from the 
same regions in the galaxy). The scatter of a factor of 2.5 is 
more than sufficient to explain the dispersion in the strength of 
the Ha emission at a fixed far-infrared luminosity. We conclude 
therefore that the large scatter in SFR(FIR)/SFR(Ha) probably 
arises from two main factors: (1) there is a range in the con- 
tinuum extinction in the more far-infrared luminous population 
and, (2) the morphological diversity of sub-mm-selected galax- 
ies - which includes a large fraction of interacting or merg- 
ing systems, sometimes with highly-obscured components may 
lead to a large variation in the in-slit Ha fluxes for these sys- 
tems (Chapman et al. 2003b; Smail et al. 2004). 

2.11. Metallicities 

Rest-frame optical emission lines from extra-galactic Hll re- 
gions provide an important diagnostic of the chemical evolu- 
tion of galaxies since their properties reflect the make-up of 
the (ISM). Moreover, feedback processes from short-lived mas- 
sive stars in high SFR galaxies may be responsible for the en- 
richment of the intergalactic medium (IGM) with heavy met- 
als over a wide range of redshifts (Heckman et al. 2000; Adel- 
berger et al. 2003). For star-forming and irregular galaxies the 
correlation between metallicity (traced by Oxygen abundances, 
O/H) and rest-frame luminosity is well established in the local 
universe and spans a factor of over 100 in (O/H) and at least 
eight magnitudes in M v . The correlation is in the sense that 
more massive (and luminous) galaxies exhibit a higher degree 
of metal enrichment (Garnett et al. 1987; Zaritsky et al. 1994; 
Garnett 2002; Lilly et al. 2003; Lamareille et al. 2004). If such 
a correlation exists at high redshift then it may be more im- 
portant for understanding the present distribution of metals in 
the universe, since the bulk of the star-formation activity in the 
most massive and active galaxies is thought to occur at z > 1 
(Blain et al. 1999b). 

Since we lack the observations necessary to calculate the 
/?23 index (Zaritsky et al. 1994) for our sample, we turn to 
the N2 index described by Storchi-Bergmann et al. (1994) 
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Fig. 8. — Comparison of the far-infrared vs. narrow-line Ha luminosities and star formation rates in our data with those in local samples. We 
include in the plot the SFRs for the SMGs/OFRGs when the broad-line Ha flux is included in the estimate (the galaxies classified as intermediate 
in §2.4 are included with the AGN classification). This is more comparable to the quantities calculated for local samples. We also include the 
previously published SMGs (Kneib et al. 2004; Dey et al. 1999; Smail et al. 2003; Frayer et al. 2003; Ivison et al. 1998) as red points ordered 
from lowest to highest Ha luminosity respectively. The dotted line represents identical estimated SFR from Ha and far-infrared. The sample is 
compared to the results from ISOCAM, including the z = 0.2-1.5 galaxies in the Hubble Deep Field-South by Franceschini et al. (2003) and the 
local luminous infra-red galaxy sample (z = 0.2-0.7) by Flores et al. (2004). We also compare the data to the Very Luminous and Ultra-Luminous 
infrared galaxy sample of local IRAS greater than 2Jy sources by Dopita et al. (2002) and Poggianti & Wu (2000). We include a representative 
error bar from the comparison samples on the left side of the plot. 



and discussed more recently by Denicolo et al. (2002) and 
Pettini & Pagel (2004). This indicator is defined at N2 = 
/og([Nll]A6584/Ha) and is calibrated to the oxygen abundance 
(O/H) via 12 + /og(0/H) = 8.90 + 0.57 x N2 (Pettini & Pagel 
2004). However, this calibration remains uncertain, and there 
are obvious drawbacks to using this index for young galaxies 
in which AGN may play a role in defining the emission-line 
characteristics and in which primary and secondary sources of 
Nitrogen production may lead to mismatches in N/O. For this 
reason we have chosen to compare populations directly using 
their N2 measurements (as opposed to (O/H) measured from a 
range of indicators) to minimise systematic uncertainties aris- 
ing from the calibration from N2 to (O/H). We therefore deter- 



mine the N2 index for all of the galaxies in our sample and give 
upper limits where the [Nil] emission line is not detected with 
sufficient significance to measure accurately. 

In Fig. 9 we show the distribution of My versus N2 for 
SMGs/OFRGs as a proxy for their Luminosity-Metallicity re- 
lation and compare this with similar observations of both local 
and high redshift galaxy populations. To calculate the observed 
rest-frame V-band luminosities of the SMGs we use the pho- 
tometry from Smail et al. (2004) and fit SED's based on their 
observed broad-band UK magnitudes, which, at z ~ 2.4, cor- 
responds to rest-frame UBR, bracketing the V-band. As Smail 
et al. (2004) show, the competing effects of dust reddening and 
fading of the young stellar populations almost cancel each other 
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FIG. 9. — Left: Metallicity-Luminosity relation for SMGs/OFRGs. Right: Line width - metallicity relation for the SMGs/OFRGs. The abundance is derived using 
the N2 index calibration given by Pettini & Pagel (2004). We show two low-redshift samples for comparison, one from the SDSS DR2 (Abazajian et al. 2004) and 
one for local luminous infrared galaxies (Annus et al. 1989). For a higher redshift comparison sample, we show the metallicities for the z ~ 2 UV-selected galaxies 
from Shapley et al. (2004) and for a near-infrared-selected sample of luminous, z. ~ 2 galaxies from van Dokkum et al. (2004). The horizontal dashed line represents 
solar metallicity. The SMGs/OFRGs exhibit a large range in My and FWHM lcs i, and many have slightly sub-solar or solar abundances. Note that these are the 
observed V-band luminosities (uncorrected for reddening; c.f. Fig. 7). 



(AdustV ~ 3-4, Af a deV ~ 4 based on fitting the broadband col- 
ors of this population). Thus, we expect that the present-day 
descendants of these galaxies would have absolute luminosities 
not very dissimilar to those we have estimated (see also Fig. 7). 
For the comparison samples we convert the rest-frame B-band 
magnitudes for the UV-selected z ~ 2 sample from Shapley et 
al. (2004) using their (R-K) colors to predict the rest-frame 
V-band magnitudes. We also convert the median, 10 and 90 
percentile trend lines from a volume-limited (z < 0.1) local 
emission-line galaxy sample from the SDSS DR2 (Abazajian et 
al. 2004; although we stress that the Sloan aperture only sam- 
ples the central ~ 3 kpc of a galaxy at z = 0.05 whereas the 
NIRSPEC slit will sample ~ lOkpc at z = 2.3). We see at least 
an order-of-magnitude range in NT for the SMGs/OFRGs (dis- 
counting obvious AGN), with little correlation between the line 
index and rest-frame luminosity. The implied median metal- 
licity is slightly below solar, and appears to be similar to that 
inferred for bright, UV-selected galaxies at z ~ 2 (Shapley et al. 
2004). 

To address the competing and uncertain effects of dust red- 
dening and passive fading of the stellar populations on the My 
versus [Nll]/Ha plot, we have also constructed a FWHMh q 
versus [Nll]/Ha diagram, adopting the FWHMh q as a crude 
proxy for the dynamical mass of the galaxies. As can be seen, 
the large dynamical masses we inferred for the SMG/OFRG 
population suggest that their present-day descendants are likely 
to be luminous and metal-rich (super-solar) systems (as shown 
by the trends seen in the SDSS dataset plotted in Fig. 9). The 
apparently modest metallicities we measure would then indi- 
cate that these systems are seen during an early phase of en- 
richment - suggesting that they are relatively youthful galax- 
ies and arguing against them undergoing a cycle of repeated 
short (~ lOMyr) bursts of star-formation over a relatively ex- 
tended period (> 1 Gyr; Smail et al. 2003). However, there are 
other possible explanations for the apparently low metallicities 
of these galaxies. First, we note that star-formation can usually 



only act to increase a galaxy's metallicity (and therefore NT), 
and in particular starburst-driven feedback mechanisms are un- 
likely to preferentially expel large quantities of metals without 
entraining and expelling associated gas. Only an infall of unen- 
riched material into the galaxy could cause the metallicity and 
effective heavy-element yield to decrease (e.g., Garnett 2002). 
Alternatively, if these systems are very young and the halos of 
the SMGs/OFRGs have yet to coalesce (as suggested by the 
clear merger/interacting morphologies of many systems), then 
the metallicities we measure may reflect those of the progenitor 
components. The similarity of the NT estimates with those for 
the UV-selected population at this epoch could then be inter- 
preted as indicating that the SMGs/OFRGs arise from mergers 
among the UV population. The slightly more evolved descen- 
dants of these mergers will be able to retain their enriched gas 
and so produce a super-solar stellar population. We stress that 
local calibrations of the NT index from Hll regions can show 
super-solar metallicities which may not reflect the global abun- 
dance in a galaxy. However, since our spectroscopic slit covers 
most of the area of our galaxies and the measured Ha fluxes 
suggest substantial star-formation, we suggest that the NT. in- 
dex should provide a fair estimate of the metallicity of the gas in 
the galaxy as a whole. We caution that the nitrogen abundances 
that we measure come from primary metal enrichment, whereas 
the local N/O versus O/H relations suggest that for metal-rich 
systems the nitrogen is secondary in origin (produced from in- 
termediate mass stars far removed from the first generation). 
The evolution of the nitrogen metallicity-luminosity relation 
might be therefore quite different from that shown for oxygen or 
other a-elements. We also note that the NT index is not an ideal 
metallicity indicator given the presence of an AGN in many of 
these galaxies, and we look forward to future studies based on 
[On], [Olll] and H/3 emission lines. 



2.12. X-ray comparisons 
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Using deep Chandra observations of the Hubble Deep Field 
(HDF), ELAIS N2 field, and SSA13 field (Alexander et al. 
2003b; Manners et al. 2003; Mushotzky et al. 2000), it 
is possible to compare the X-ray and Ha properties of the 
SMGs/OFRGs which overlap with the Chandra coverage. Of 
the 18 sources which were covered by Chandra, nine were de- 
tected in the 2-8 keV hard X-ray band, mostly in the HDF ow- 
ing to the much deeper X-ray observations available for that 
field. We convert the observed 2-8 keV flux to a rest-frame 
2-10keV luminosity using Lx = 4-irdlfx(l +z) r ~ 2 , which takes 
into account the ^-correction (Alexander et al. 2003b), assum- 
ing a spectral index T = 2. 

In Fig 10 we compare the Ha luminosities from the 
SMGs/OFRGs to their X-ray luminosities and contrast these 
with a local sample of /RAS-selected Seyfert galaxies anal- 
ysed by Ward et al. (1988). We also overlay two lines show- 
ing the correlation of L(2-ioKeV) versus L(FIR) from Ranalli et 
al. (2003). The first line simply converts their L(FIR) to L(Ha) 
assuming that the relation in §2 holds, and the second line as- 
sumes a further 10 times suppression of the Ha emission (rela- 
tive to the far-infrared) as indicated by Fig. 8. 

The distribution of the SMGs/OFRGs in Fig. 10 roughly fol- 
lows that seen for the Ward et al. (1988) sample - irrespective 
of the spectral classification of the SMGs/OFRGs. However, 
we note that the correlation for local star-forming galaxies from 
Ranalli et al. (2003) - if scaled for the relative underluminosity 
of Ha in the SMGs/OFRGs - can explain the properties of the 
least X-ray luminous galaxies, these comprise roughly half of 
our sample. 
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FIG . 10.- Distribution of Ha and X-ray luminosities for the SMGs/OFRGs in 
our sample as compared to the local /SAS-selected Seyfert galaxies from Ward 
et al. (1988). For the SMGs/OFRGs sample we plot the total Ha luminosity 
from both narrow and broad components (if present). We also compare the data 
to a sample of nearby star-forming galaxies by Ranalli et al. (2003), shown by 
a dotted line (converting their L(FIR) to L(Hcr) assuming the correlation in §2 
holds) and include the same line but assuming a 10 times suppression of the Ha 
emission (relative to the far-infrared) as indicated in Fig. 8, shown by a dashed 
line. The plot shows that some SMGs/OFRGs which are identified as star- 
forming in the near- infrared (with FWHM < SOOkms -1 and [Nll]/Ha S, 0.7) 



show 2-10 keV X-ray fluxes consistent with AGN luminosities. It is also in- 
teresting to note that some of the SMGs/OFRGs which are spectroscopically- 
classified as AGN from their near-infrared spectra are not detected in the hard 
X-ray band. 

The most intriguing galaxies in Fig. 10 are those which have 
high X-ray luminosities (logL A > 36.5) but have restframe op- 
tical spectral classifications of star-formation - the X-ray lumi- 
nosity in these galaxies probably arises from an AGN which is 
so highly obscured that either it is hidden in the rest-frame opti- 
cal spectra (given their modest signal to noise) or we missed the 
AGN component with our spectroscopic slit in these systems. 
The presence of a broad-line component in the composite star- 
forming spectrum in §3.2 would support the former interpreta- 
tion. 

We also point out that there are six SMGs/OFRGs which are 
not detected to a flux limit of 1.4 x 10" 19 Wm" 2 in the 2-Ms 
Chandra observations of the HDF (two) or 2.2 x 10" I8 Wm" 2 
in the ELAIS N2 (four) field. While three of these are clas- 
sified as SB from their restframe optical spectra, three others 
show some signs of AGN in the near-infrared spectra (with 
broad Ha or large [Nll]/Ha ratios). The X-ray limits on two 
of these (SMM J163650.43 and SMM J163706.51) are partic- 
ularly stringent - yet they show clear AGN signatures in their 
near-infrared spectra. The detection of broad Ha and lack of X- 
ray emission are difficult to understand - with these galaxies be- 
ing almost an order of magnitude fainter in the X-ray waveband 
than expected from their Ha luminosities. Given the amount of 
gas available to fuel an AGN in these systems, we conclude that 
they are likely to be intrinsically low luminosity and probably 
have low mass central black holes. 

3. CONCLUSIONS 

We present the results of near-infrared spectroscopic and 
narrow-band detections or limits on the Ha emission from a 
sample of 30 ultraluminous, dusty galaxies at z ~ 1.4-2.7. The 
majority of these galaxies come from sub-mm/mm-surveys, 
with a small number identified as probable hot, luminous far- 
infrared sources from their radio emission. We see no differ- 
ence between the properties of these two samples in any of our 
diagnostic diagrams - supporting the claimed similarity of the 
two populations (Chapman et al. 2004b). 

We identify the Ha emission in the near-infrared spectra and 
use the spectra to classify AGN by flagging galaxies with large 
[Nll]/Ha ratios (> 0.7) and/or large Ha FWHM (> SOOkms" 1 ). 
We find that the ratio of AGN- to star formation- dominated 
galaxies from the rest-frame optical spectroscopy is roughly 
40%:60% - similar to the proportions estimated from their UV 
spectra by Chapman et al. (2003a, 2003b, 2004c). By con- 
structing a rest-frame composite spectrum for the entire sam- 
ple, we find that average [Nll]/Ha ratio is 0.42 ± 0.05, which 
suggests that the composite is star-formation- rather than AGN- 
dominated, although it has an underlying broad-line compo- 
nent. Furthermore the composite spectra show both NaD ab- 
sorption and [Sll] emission features, although we find no ev- 
idence of velocity offsets between these features and the Ha 
emission. The [Sll]/Ha emission ratio in the composite spec- 
trum is 0.10 ± 0.04 - indicating that the spectral properties of 
our sample are comparable to a LINER or Hll region - sim- 
ilar results are found from the spectral classifications of local 
ULIRGs (Veilleux et al. 1987, 1995). 

We also derive the composite spectrum for those galaxies 
which, individually, show no signs of an AGN in their near- 
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infrared spectra. This composite has an average [Nll]/Ha flux 
ratio of 0.19 ± 0.05 and FWHM rest of 400 ± 70^8"'. How- 
ever, the most striking result is that the composite spectrum ap- 
pears to show an underlying broad Ha line with a broad/narrow 
Ha flux ratio of 0.45 ±0.20, suggesting that even these galaxies 
may host a low-luminosity AGN which is undetectable in our 
modest signal-to-noise spectroscopy. 

In seven of the systems with spectroscopic observations we 
find velocity structure in the Ha emission line. By com- 
parison with high-resolution broadband imaging, we identify 
these galaxies as multi-component (probably interacting) sys- 
tems with typical velocity offsets between components of 100- 
600 km s" 1 . This is not surprising since locally many far- 
infrared luminous galaxies appear to be disturbed/interacting 
systems. Assuming that these are merging systems with ran- 
dom orientations of their orbits on the sky, we estimate a typical 
mass of 1 .5 ± 0.9 x 1O U M . We obtain a similar estimate from 
the Ha line widths of the whole sample. These estimates are 
comparable to the dynamical mass estimates from CO observa- 
tions of a subset of these systems. 

In all of the galaxies, we have attempted to deconvolve any 
broad component to the Ha line (which comes from an accre- 
tion disk around a central super-massive black hole) from the 
narrow-line Ha flux (which comes from the star-forming re- 
gions). Using the narrow-line Ha flux, we compare the SFRs 
of the SMGs/OFRGs as compared to the SFR derived from 
the far-infrared emission. The SFR(Ha) versus SFR(FIR) cor- 
relation shows a large scatter, with the SFR(Ha) typically a 
factor of ten less than what we would expect from their far- 
infrared luminosities. Most of this scatter, however, can be 
explained by the reddening in these systems (estimated from 
their broad-band photometry in Smail et al. (2004)). The sup- 
pression of the Ha flux is therefore attributed to both heavily 
obscured galaxies and a diverse range of morphologies. The 
average SFR derived from Ha for the SMGs/OFRGs in our 
sample is 94 ± 20M Q yr _1 (uncorrected for extinction). Since 
the continuum extinction correction at 6563A is ^2.5 magni- 
tudes, the total inferred SFR of these galaxies is expected to 
be ^1000M Q yr _1 - comparable to that seen in the far-infrared. 
We also compare the SFR properties for our high-redshift sam- 
ple to local Very/Ultra-Luminous galaxy samples and find that 
the scatter within our high redshift SMGs/OFRGs and the dis- 
tribution of equivalent widths are comparable to those of local 
Ultra/Very luminous galaxy samples. This suggests a similar 
range of obscured/unobscured activity in the distant population 
to that seen locally, although with a higher proportion of the 
star formation completely obscured from view. 

Using the N2 index, we have investigated the chemical abun- 
dances of these galaxies and find that the N2 indices for the 
SMGs/OFRGs suggest that they have slightly sub-solar metal- 
licities, similar to recent results from UV selected galaxies at 
these early epochs (Shapley et al. 2004). However, we note that 
the N2 indicator may not be a reliable metallicity indicator for 
this population - especially in the presence of an AGN. We find 
that the SMGs/OFRGs in our sample display a large range in 
[Nll]/Ha versus My or FWHM rest (which we use as a proxy for 
their dynamical masses). 

For the galaxies which are in the HDF, ELAIS N2 field and 
SSA13 field we use their X-ray properties in order to further 
classify the galaxies. The SMGs/OFRGs classified as AGN 
on the basis of their near-infrared spectra broadly follow the 
correlation seen between X-ray and Ha luminosities for local 



Seyfert 2 galaxies (Ward et al. 1988). We find that a subset of 
galaxies which are spectroscopically-classified as star-forming 
in the near-infrared have high X-ray luminosities, suggesting 
that they host highly obscured AGN. Likewise, a small num- 
ber of near-infrared spectroscopically-classified AGN are un- 
detected in deep Chandra observations. We conclude that these 
galaxies are likely to be intrinsically low luminosity and prob- 
ably have low-mass central black holes. 

With observations of the Ha emission from these 
SMGs/OFRGs at z ~ 2 we can at last start to directly compare 
these galaxies to similarly distant UV-selected systems. We find 
that as expected the SFR(Ha) for our sample is nearly an or- 
der of magnitude higher than that found in z ~ 2 UV-selected 
galaxies by Erb et al. (2003) (who find an average SFR(Ha) 
of 21 ± 3M Q yr _1 - even after their sample is corrected for ex- 
tinction). Similarly, the line widths and dynamical information 
suggest that the halos of a typical SMG/OFRGs may be upto 
5 times more massive than the UV-selected population. We 
also find a higher rate of AGN activity in the SMG/OFRGs - 
suggesting the presence of actively fueled and growing super- 
massive black holes in these galaxies. However, somewhat 
surprisingly, we find similar metallicities for the UV-selected 
and more massive far-infrared luminous populations. We sug- 
gest that this may be explained if the SMG/OFRGs are rel- 
atively youthful, with their deepening potential wells not yet 
sufficiently organised to retain a larger fraction of the enriched 
material from their star formation activity. 

Overall our observations suggest that the high redshift 
SMG/OFRG population shares many of the characteristics 
of similar (but somewhat less luminous) far-infrared galaxies 
identified in the local Universe. This includes the Ha equiv- 
alent widths, the proportion of obvious AGN and the typical 
optical spectral classification. Previous work has demonstrated 
the preponderance of merger-like morphologies in the two pop- 
ulations and the similarity of their restframe optical luminosi- 
ties. Yet there remain differences, with proportionally more 
highly-obscured activity in the high-redshift population, appar- 
ently larger dynamical mass, lower metallicities and higher gas 
fractions on 10-kpc scales. 

We conclude that the SMG/OFRGs in our sample represent 
a population of young, massive merging/interacting systems, 
the results of which cause high instantaneous bursts of (highly 
obscured) star-formation and actively fueled AGN activity. Al- 
though these bursts are brief, they can form all of the stars in 
an L* galaxy and in doing so will raise the metallicity of these 
systems closer to that required by observations of their likely 
present-day descendants: luminous elliptical galaxies. 
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